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Indications for Criticality at Zero Curvature in a 4d Regge Model of 
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We re-examine the approach to four-dimensional Euclidean quantum gravity based on the Regge calculus. A 
cut-off on the link lengths is introduced and consequently the gravitational coupling and the cosmological constant 
become independent parameters. We determine the zero curvature, (R) = 0, line in the coupling constant plane 
by numerical simulations. When crossing this line we find a strong, probably first order, phase transition line 
with indications of a second order endpoint. Beyond the endpoint the transition through the (R) = line appears 
to be a crossover. Previous investigations, using the Regge or the Dynamical Triangulation approach, dealt with 
a limit in which the first order transition prevails. 



1. Introduction 

Euclidean quantum gravity in four dimensions 
is formally defined by the path integral 
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The Einstein-Hilbert action Seh\9\ IS a function 
of the 4-geometry g, SehM = -Lp 2 R[g]+AV[g}. 
The Planck length Lp determines the gravita- 
tional coupling constant for the total curvature 
R[g] and the total 4- volume V[g] enters the ac- 
tion if the cosmological constant A is non-zero. 

In order to define the integral (1), we consider 
a lattice approximation, so that the path-integral 
becomes the limit of a series of well-defined and 
finite integrals. A direct route for such an ap- 
proximation is provided by the Regge calculus [1] , 
which is suitable for simulations [2,3]. One con- 
siders a simplicial lattice of fixed connectivity, 
which subdivides a 4-manifold of given topology 
in N4 4-simpliccs, A3 tetracder, N2 triangles, Ni 
links and No vertices. As we assign a length 
to each link, a 4-geometry is defined in the Eu- 
clidean sector, if the triangle inequalities are ful- 
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filled for every (sub-)simplex. In the following we 
will use the squared lengths qi associated with 
each link I as the basic variables and a given set 
of such variables defines a 4-geometry [q] for a 
given simplicial lattice. Variation of the vari- 
ables qi allows us to approximate the path in- 
tegral (1). One can calculate the Regge action 
R[q] = J2t as a sum over au triangles t with 
area A t multiplied by their associated deficit an- 
gles St. The 4- volume is calculated as a sum over 
the volume of all 4-simplices s, V[q] = ^2 S V S . 
We then approximate the gravitational action by 
S eh -> S[q] = -L p 2 R[q] + AV[q]. It has been 
demonstrated how the geometrical variables of 
the Regge calculus converge towards the corre- 
sponding variables of General Relativity in the 
classical continuum limit [1,4]. It requires an up- 
per bound on the link lengths, qi < Q, and a 
limit F on the 'fatness' <j) s =V S / max(g ; 2 ) of each 
4-simplex, <p s > F. The continuum limit is then 
reached for Ao — > 00, Q — > with F > (Cheeger 
et al. [4]). 

In order to replace the path integral (1) by a 
finite integral on the simplicial lattice, 



Dqe 



S[q] 



(2) 



2 



we need to define the integration measure. A sim- 
ple choice is the uniform measure 

Dq= [ T\d qi e F [q] (3) 

Jo , 

where the function <d F equals 1 when all triangle 
inequalities arc fulfilled and <f> s > F for all 4- 
simplices; it is otherwise. The integration limit 
Q is necessary for a well-defined integral, but it 
has no physical meaning in itself. Indeed, simple 
re-scaling qi — > Qqi formally removes it and leads 
to the integral 

Z(P,\)= ! T\d qi Q F [q}e 0R ^- xv ^ (4) 
Jo , 

where (3 = QL p 2 and A = Q 2 A. 

In previous numerical investigations [5,6] the 
upper limit Q was infinite which corresponds to 
the limit [3 — > oo and A — > oo in the integral (4). 

In the following we will numerically evaluate 
the integral (4) in order to determine the phase 
structure by varying the coupling parameters (3 
and A. Our results suggest a physically interest- 
ing continuum limit for finite values of the cou- 
pling parameters. 

2. Results 

The numerical computations have been per- 
formed for the integral 

Z(P,\)= f " l[d qi e F [q}e^-^ (5) 
Jo l 

with c = 10 instead of (4) with c = 1. This means 
that the coupling constants and expectation val- 
ues need to be adjusted accordingly, e.g. /3 — ► c/3, 
A — ► c 2 A etc. We set F = 10~ 6 but our results 
seem to be independent of the specifc value for 
a relatively wide range F = 10~ 8 to 10~ 5 . We 
used a hypercubic regular triangulation of the 4- 
torus [2,3], with N a = 4 4 ,4 3 x 10, 6 4 and 8 4 ver- 
tices. 

From simulations with A = 1,2,3,4 and vary- 
ing f3 we find expectation values < R >, which 
are (almost) zero for /?o(A) values; /3 initially de- 
creases for increasing A. Then, as depicted in fig- 
ure 1, we find that /3 increases for A = 4, 5, 6, 7. 
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Figure 1. Hysteresis loops. 
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Figure 2. Function /3o(A) for which < R >= 
holds. 

Furthermore, the transition from < R >< to 
< R >> is no longer defined clearly, because 
we have to deal with a significant hysteresis for 
large A. The line (3o(X) is depicted in figure 2 
for different lattice size. The curve /3o(A) is in- 
teresting, since a vanishing curvature is a desired 
ingredient for a physical continuum limit. 

As we approach the turning point A c ~ 4.0 from 
smaller A, the suszeptibility \ q = N^ 1 (< q 2 > 
— < q > 2 ), with q = N^ 1 J2i1h increases signifi- 
cantly; this increase is more pronounced for larger 
lattice size as depicted in figure 3. 

Next, we performed calculations varying A at 
= 0.175 and 0.172 for the lattice size 4 3 x 10. 
As expected, we find that the curvature < R > in- 
creases initially, exhibiting a transition from neg- 
ative to positive values, but when A increases fur- 
ther a second transition back to negative values is 
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Figure 3. Susceptibility Xq as function of A at (3q 
for increasing lattice size. 
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Figure 4. Susceptibility xr as function of A for 
two (3 values. 

seen. We depict xr = ^2 _1 (< R 2 > - <R > 2 ), 
measuring curvature fluctuations, in figure 4; for 
(3 = 0.175 one can distinguish two peaks cor- 
responding to the two transitions. These peaks 
merge into one at (3 = 0.172. 

We also measured the 2-point correlation func- 
tion K(d) =< q qd > — < qo > 2 along the elon- 
gated direction on the 4 3 x 10 lattice for /3o(A) 
and found that the function K{d) is significantly 
different from zero for d > 1. Negative values 
alternate with positive values, which indicates a 
non-trivial interaction. The correlations increase 
for A — ► A c . 

Although we are not able to provide a com- 
plete finite size scaling analysis, our data suggest 
the following: In the region A > A c , with A ~ 4.0, 
one deals with a strong phase transition, which is 



likely of first order. In the region A < A c we are 
currently unable to distinguish between a cross- 
over transition or a higher-order phase transition. 
However the end-point A c itself is a natural can- 
didate for a second-order phase transition. The 
increase of susceptiblity Xq an( f the link correla- 
tions K(d) approaching A c from smaller values 
confirm this picture. 

3. Summary and Conclusions 

Our simulations differ from those of the pre- 
vious literature by introducing the limit Q in 
the functional integral (3). While a simple re- 
scaling removes Q formally from the partition 
function (4), our model remains more general 
than those simulated in the literature, which are 
recovered in our limit of large A. Remarkably, the 
new model appears to have a second order fixed 
point in the A — (3 plane, which deserves further 
study. 
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